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4 Optimum Reception in Additive White Gaussian 
Noise (AWGN) 

In this chapter, we derive the optimum receiver structures for the modu- 
lation schemes introduced in Chapter 3 and analyze their performance. 

4.1 Optimum Receivers for Signals Corrupted by AWGN 

■ Problem Formulation 

— We first consider memoryless linear modulation formats. In 
symbol interval 0 < t < T, information is transmitted using 
one of M possible waveforms s m (£), 1 < m < M. 

— The received passband signal r{t) is corrupted by real-valued 
AWGN n(t): 

r(t ) = s m (t) + n(t), 0 < t < T. 


Sm(t) 


n(t), N 0 / 2 




r(t) 


— The AWGN, n(£), has power spectral density 

rwi 


^nn(J) = ~y 


Hz 


— At the receiver, we observe r(t) and the question we ask is: 

What is the best decision rule for determining s m (t)l 

— This problem can be equivalently formulated in the complex 
baseband. The received baseband signal r&(£) is 

Tb(t) Sbm(t) T zffi) 


where z(t) is complex AWGN, whose real and imaginary parts 
are independent. z(t) has a power spectral density of 


$zz(f) = 


"W" 

Hz 


s bm(t ) 


z(t), No 


-e- 


nit) 


Strategy: We divide the problem into two parts: 

1. First we transform the received continuous-time signal r(t) (or 
equivalently r&(£)) into an TV-dimensional vector 

r = [n r 2 ... r N ] T 

(or r&), which forms a sufficient statistic for the detection of 
s m (t) (s^ m (t)). This transformation is referred to as demodu- 
lation. 

2. Subsequently, we determine an estimate for s m (t) (or s& m (£)) 
based on vector r (or r &) . This process is referred to as detec- 
tion. 
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r(t) 


Demodulator 


r 


Detector 


Decision 


4.1.1 Demodulation 

The demodulator extracts the information required for optimal detec- 
tion of s m (t) and eliminates those parts of the received signal r(t) that 
are irrelevant for the detection process. 

4. 1.1.1 Correlation Demodulation 

■ Recall that the transmit waveforms (s m (£)} can be represented by 
a set of N orthogonal basis functions /&(£), 1 < k < N. 

■ For a complete representation of the noise n(t), 0 < t < T, an 
infinite number of basis functions are required. But fortunately, 
only the noise components that he in the signal space spanned by 

1 < k < N, are relevant for detection of s m (t). 

u We obtain vector r by correlating r(t) with 1 < k < N 


T 


T 



0 


0 


T 


T 



Smk T ^ ki 


1 < k < N 



■ r(t) can be represented by 

N N 

r(t) = mkfk(t ) T ^ ^ 'H'kfkjt) + Tl (t) 
k= 1 k= 1 

N 

= ^2nf k (t) + n\t), 
k—1 

where noise n'(t) is given by 

N 

n'(t ) = n(t) - y^rifc/fc^) 

jfe=i 

Since n'(t) does not lie in the signal space spanned by the basis 
functions of s m (£), it is irrelevant for detection of s m (t). There- 
fore, without loss of optimality, we can estimate the transmitted 
waveform s m (t) from r instead of r(t). 


Properties of n\~ 

— rik is a Gaussian random variable (RV), since n(t) is Gaussian. 

— Mean: 


£{n k } m £<^J n(t)ft(t)dt 

T 

= j S{n{t)}f* k (t) At 


0 

= 0 


— Covariance: 


£{«*<>} = £ j yj n{t)f k {t)dtj yj n(i)/*(i) 

T T 

= J J f{n(t)n*{T)\ fi(t)f m (T) At At 

0 0 !$S(t-T) 

T 

= Y J 


d t 


No n 

= —S[k-m] 


where 5[k] denotes the Kronecker function 
1, k = 0 
0, 0 

We conclude that the N noise components are zero-mean, mu- 
tually uncorrelated Gaussian RVs. 


S[k] = 


■ Conditional pdf of r 

r can be expressed as 


r = s m + n 


with 


Sm [^ml Sm2 ■ • • ^miv] 5 

n = [m n 2 ... n N ] T . 

Therefore, conditioned on s m vector r is Gaussian distributed and 
we obtain 

p{r\s m ) = p n (r — s m ) 

N 

= 1 J. Pnfrk ^ mk)i 

k= 1 

where p n (n) and p n {^k) denote the pdfs of the Gaussian noise 
vector n and the components n k of n, respectively. Pni^k) is 
given by 


PniP'k ) — 


1 


a/ttTVo 


exp 1 -N 0 


since n k is a real-valued Gaussian RV with variance crj; = -j-. 
Therefore, p(r\s m ) can be expressed as 

N / /_ n 2 > 


P{r\s m ) = \ 


1 


/¥Nr 


exp 


(tyNq) N /‘ 


: exp 


{j*k Smk) 

A^o 

/ * \ 

/ ^ (j'k — Smk) 

k= 1 

Nn 


1 < m < M 


/ 


p(r\s m ) will be used later to find the optimum estimate for s m (or 
equivalently s m (t)). 

■ Role of n'(t) 

We consider the correlation between r & and n'(t)\ 

£{n'(t)r* k } = £{n'(t)(s mk + n k )*} 



= 0 


We observe that r and n'(t) are uncorrelated. Since r and n'(t) 
are Gaussian distributed, they are also statistically independent. 
Therefore, n'(t) cannot provide any useful information that is rele- 
vant for the decision, and consequently, r forms a sufficient statis- 
tic for detection of s m (t). 
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4.1. 1.2 Matched-Filter Demodulation 

■ Instead of generating the {r^} using a bank of N correlators, we 
may use N linear filters instead. 

■ We define the N filter impulse responses hk(t) as 

fk( T -a o<t<r 

where /&, 1 < k < N, are the N basis functions. 

■ The output of filter hk(t) with input r(t) is 

t 

Vk(t) = J r(r)hk(t — t) dr 

o 

t 

= J r{r)fl(T -t + r)dr 
0 

■ By sampling yk(t) at time t = T, we obtain 

T 

Vk (T) = J r(r)/ t *(r) dr 
0 

= rfc, 1 < k < N 


This means the sampled output of hk(t) is r^. 



t = T 


General Properties of Matched Filters MFs 

— In general, we call a filter of the form 

h(t) = s*(T — t) 

a matched filter for s(t). 

— The output 

t 

y(t) = J s(t)s*(T — t + t) dr 
o 

is the time-shifted tim^autocorrelation of s(t). 
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Example: 



- MFs Maximize the SNR 

* Consider the signal 

r(t) = s(t) + n(t), 0 < t < T, 

where s(t) is some known signal with energy 
T 

E = J |s(t )| 2 d t 

o 

and n(t ) is AWGN with power spectral density 
^nn(I) = -j- 
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* Problem: Which filter h(t) maximizes the SNR of 

y(T) = h(t ) * r(t) 

t—T 

* Answer: The matched filter h(t) = s*(T — t)\ 

Proof. The filter output sampled at time t = T is given by 

T 

y(T ) = J r(r)h(T — r) dr 
o 

T T 

= J s(r)h(T — t) dr + J n(r)h(T — r) dr 


ys( T ) vn( t ) 

Now, the SNR at the filter output can be defined as 

SNR = Jy^_ 

n\y N m 

The noise power in the denominator can be calculated as 
£ {12/jvCOI 2 } = £ |^y~ n(r)h(T - r) dr j n(r)h(T - t) dr ^ 

T T 

= J J Z {n(r)n* (t)} h(T — r)h* (T — t) dr < 

o 0 ^S(r-t) 

= Y J \h(T-r)\ 2 dr 


-At 


Therefore, the SNR can be expressed as 


SNR = 


f s(r)h(T — r) dr 


f /KT-T)) 2 dT 


From the Cauchy-Schwartz inequality we know 

J s(r)h(T — r) dr < J |s(r)| 2 dr- J \h(T — r)| 2 dr, 

0 0 0 
where equality holds if and only if 

h(t) = Cs*(T — t). 

( C is an arbitrary non-zero constant). Therefore, the max- 
imum output SNR is 

2 


SNR = 


/|s(r)| 2 dr 


fi/|s(r)| 2 dr 


2 

N o 


1 

J kWI 5 


dr 


-4 

which is achieved by the MF h(t) = s*(T — t). 


- Frequency Domain Interpretation 

The frequency response of the MF is given by 


H(f) = F{h(t)} 

OO 

= J s*{T-t)e~ j27,ft dt 

— OO 
OO 

= J s*{T)e j27rfT e- j2 * fT dr 

— OO 



= e~ j 2 nfT S*(f) 


Observe that H(f) has the same magnitude as S(f) 

\H(f)\ = \S(f)\. 

The factor e~^ 2n ^ T in the frequency response accounts for the 
time shift of s*(—t) by T. 
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4.1.2 Optimal Detection 

Problem Formulation: 

■ The output r of the demodulator forms a sufficient statistic for 
detection of s m (t) (s m ). 

■ We consider linear modulation formats without memory. 

■ What is the optimal decision rule ? 

u Optimality criterion: Probability for correct detection shall be 
maximized, he., probability of error shall be minimized. 

Solution: 

■ The probability of error is minimized if we choose that s ^ which 
maximizes the posteriori probability 

P(srh\r), in = 1, 2, . . . , M, 

where the ” tilde” indicates that s ^ is not the transmitted symbol 
but a trial symbol. 
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Maximum a Posteriori (MAP) Decision Rule 

The resulting decision rule can be formulated as 

m = argmax {P(s^|r)} 


where rh denotes the estimated signal number. The above decision 
rule is called maximum a posteriori (MAP) decision rule. 

■ Simplifications 

Using Bayes rule, we can rewrite P(s^|r) as 

p(r) 

with 


p{ r \ s fh)- Conditional pdf of observed vector r given s^- 

— P(s ? f i ): A priori probability of transmitted symbols. Nor- 
mally, we have 

p (*rn) = jj, 1 < m < M, 

i.e., all signals of the set are transmitted with equal probability. 

— p{r)\ Probability density function of vector r 

M 

p{r ) = ^ / p{r\s m )P(s m ). 

m = 1 

Since p(r) is obviously independent of we can simplify the 
MAP decision rule to 

m = argmax (p(rfy^)P(s^)} 
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Maximum-Likelihood (ML) Decision Rule 

■ The MAP rule requires knowledge of both p(r\8jh) and P(s^). 

■ In some applications P(sm) is unknown at the receiver. 

■ If we neglect the influence of P(s^), we get the ML decision rule 


m = argmax {p(r \ 

m 


■ Note that if all s m are equally probable, he., P(s^) = 1/M, 1 < 
m < M, the MAP and the ML decision rules are identical. 

The above MAP and ML decision rules are very general. They can be 
applied to any channel as long as we are able to find an expression for 
p{r\s, h ). 
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ML Decision Rule for AWGN Channel 


■ For the AWGN channel we have 


/ JL 


^ ^ | f'k Srhk\ 


\ 



k—1 


N 0 


1 < rh < M 




■ We note that the ML decision does not change if we maximize 
\n(p(r\sfh)) instead of p(r\sm) itself, since ln(-) is a monotonic 
function. 

■ Therefore, the ML decision rule can be simplified as 

m = argmax {p(r\s r f l )} 


■ Interpretation: 

We select that vector s ^ which has the minimum Euclidean dis- 
tance 


rh 


argmax (ln(p(r|s^))} 


rh 



argmin {|| r - s^|| 2 } 


rh 


D(r , Sm) = \ \r - s^|| 
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from the received vector r. Therefore, we can interpret the above 
ML decision rule graphically by dividing the signal space in deci- 
sion regions. 

Example: 

4QAM 


m = 2 

• 

• 

rh — 1 



• r 


m = 3 

• 

• 

rh = 4 
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■ Alternative Representation: 

Using the expansion 

||T--s ft || 2 = ||r|| 2 -2Re{r«s ft } + ||s ft || 2 , 

we observe that ||r|| 2 is independent of s,„ . Therefore, we can 
further simplify the ML decision rule 

m = argmin {|| r - s^|| 2 } 


with 


= argmin {-2Re{r*s^}+ ||s^|| 2 } 

m 

= argmax {2Re{r • - Hs^ll 2 } 

m 


= argmax 


Re 


jr{t)sUt) 


d t 



T 

Em = J |s*(f)| 2 df. 

0 

If we are dealing with passband signals both r(t) and sjt (£) are 
real-valued, and we obtain 



[ 



Example: 

M - ary PAM transmission (baseband case) 

The transmitted signals are given by 

Sbm(t) = 

with A m = (2m — 1 — M)d, m = 1, 2, ... , M, 2d: distance 
between adjacent signal points. 

We assume the transmit pulse g(£) is as shown below. 

k9(t) 


a 


[ 


ml 

In the interval 0 < t < T, the transmission scheme is modeled as 



1. Demodulator 

— Energy of transmit pulse 

T 

E g = f \g(t)\ 2 dt = a 2 T 
0 


— Basis function f(t) 


m = 




9(t) 


= } v'T’ 

o, 


0 < t < T 
otherwise 


— Correlation demodulator 
T 

r b = J r b (t)f*(t)dt 
0 


= 

0 

T T 

= -j= J Sbm{t)dt + -^= J z(t)dt 


s bm 

= Sbm H - % 


s hm is given by 


Sbm 


VT 


i 

J Amd(t) 


d t 
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A m adt 


-Tr!- 

0 

= a\[r A rn = yj 'E gA m . 
On the other hand, the noise variance is 

= £{M 2 } 

T T 


j, J J £{z(t)z*{T)\ dtdT 


0 0 N 0 6(t-T ) 


= No 


J73] 


- pdf p(r b \A m )\ 


p{r b \A m ) = — exp I “ 


| n> \j EgA m | 

No 


2 . Optimum Detector 

The ML decision rule is given by 

m = argmax{ln(p(r fe |M^))} 

m 

= argmax{— |r 6 - y/EgA^ 2 } 

m 

= argmin{|r fe - y^^l 2 } 

m 

Illustration in the Signal Space 


TO = 1 


TO = 2 


TO = 3 


TO = 4 
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4.2 Performance of Optimum Receivers 

In this section, we evaluate the performance of the optimum receivers 
introduced in the previous section. We assume again memoryless mod- 
ulation. We adopt the symbol error probability (SEP) (also referred 
to as symbol error rate (SER)) and the bit error probability (BEP) 

(also referred to as bit error rate (BER)) as performance criteria. 

4.2.1 Binary Modulation 

1. Binary PAM (M = 2) 

■ From the example in the previous section we know that the 
detector input signal in this case is 

r b = ^/E~ g A m -hz. m = 1, 2, 

where the noise variance of the complex baseband noise is cr 2 z = 

N 0 . 

■ Decision Regions 


TO = 1 

• 

—y/Egd 


to = 2 


y/Egd 
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■ Assuming s\ has been transmitted, the received signal is 

n = - VEgd + 2 ; 

and a correct decision is made if 

tr < 0, 

whereas an error is made if 


tr > 0 , 

where tr = Re{r^} denotes the real part of r. tr is given by 

f‘li = + Zr 

where zr = Re{z} is real Gaussian noise with variance a 2 = 

Ao/2. 

■ Consequently, the (conditional) error probability is 


P{e\si) = J Pr R (r R \si)dr R . 
0 

Therefore, we get 



where we have used the substitution x = y/2(rn+ y/E~ g d) / v^Vq 
and the Q-function is defined as 



The BEP, which is equal to the SEP for binary modulation, is 
given by 

P b = P(s 1 )P(e\s 1 ) + P(s 2 )P(e\s 2 ) 

For the usual case, P(si) = P(s 2 ) = 5 , we get 

a = \p{e |si) + \p{e\s 2 ) 

= -P(e|«i), 

since P(e|si) = P(e|s 2 ) is true because of the symmetry of the 
signal constellation. 

In general, the BEP is expressed as a function of the received 
energy per bit E^. Here, E & is given by 

E b = £{\^E g A m \ 2 } 

- £,(!(-< I) 2 + !(<!)’) 

= E g d 2 . 

Therefore, the BEP can be expressed as 

Pb = Q 



■ Note that binary PSK (BPSK) yields the same BEP as 2PAM. 


2. Binary Orthogonal Modulation 

■ For binary orthogonal modulation, the transmitted signals can 
be represented as 

VE h - 
0 
0 

VE h _ 

■ The demodulated received signal is given by 

V ~Eb + n l 

r = 

n 2 

and 


si = 

s 2 = 


ni 

j* — 

_ y/Eb + ^2 _ 

if Si and s 2 were sent, respectively. The noise variances are 
given by 

<4 = £{«?} = °i 2 = £{ n l} = Y’ 

and ni and n 2 are mutually independent Gaussian RVs. 


178 | 

Decision Rule 

The ML decision rule is given by 

m = argmax {2r • - ||s^|| 2 } 

rh 

= argmax {r • , 

rh 

where we have used the fact that 1 1 s m|| 2 = El, is independent 
of rh. 

Error Probability 

— Let us assume that m = 1 has been transmitted. 

— From the above decision rule we conclude that an error is 
made if 


r • Si < r • S 2 

— Therefore, the conditional BEP is given by 

P(e\si) = P{r • s 2 > r • Si) 

= P{y/E h n 2 > E b + 

= P{ n 2 ~ n x > y/Et) 
x 

Note that X is a Gaussian RV with variance 

a 2 x = £ {\n 2 -ni| 2 } 

= £ {M 2 } - 2£ {nin 2 } + £ {|ni| 2 } 
= N 0 


Therefore, P(e\s{) can be calculated to 



Finally, because of the symmetry of the signal constellation 
we obtain = P(e|si) = P(e|s2) or 



Comparison of 2PAM and Binary Orthogonal Mod- 
ulation 

— PAM: 

Pb = Q 

— Orthogonal Signaling (e.g. FSK) 

Pb = Q 

— We observe that in order to achieve the same BEP the E &- 
to No ratio (SNR) has to be 3 dB higher for orthogonal sig- 
naling than for PAM. Therefore, orthogonal signaling (FSK) 
is less power efficient than antipodal signaling (PAM). 





E b /N 0 [dB] 


2PAM 


2FSK 




-• #- 

—y/Eb \fE~b 


V 

y/Ei 


We observe that the (minimum) Euclidean distance between 
signal points is given by 


for 2PAM and 2FSK, respectively. The ratio of the squared 
Euclidean distances is given by 


Since the average energy of the signal points is identical for 
both constellations, the higher power efficiency of 2PAM can 
be directly deduced from the higher minimum Euclidean 
distance of the signal points in the signal space. Note that 
the BEP for both 2PAM and 2FSK can also be expressed 



and 


= ^2E b 



as 



[ 


- Rule of Thumb: 

In general, for a given average energy of the signal points, 
the BEP of a linear modulation scheme is larger if the min- 
imum Euclidean distance of the signals in the signal space 
is smaller. 


4.2.2 M- ary PAM 

■ The transmitted signal points are given by 

s bm = * s [& g A rn , 1 <m< M 


with pulse energy E g and amplitude 

A m = (2m — 1 — M)d, 1 < m < M. 

i Average Energy of Signal Points 


E s = 


1 

771=1 

1 M 

-E g d 2 J2^m-l~Mf 

777=1 

,2 r M M 

4 y^m 2 — 4(M + 1) y^m + 


E g d 2 


M 


M 2 - 1 


777=1^ 

jAf(Af+l)(2Af+l) 

■ d 2 E„ 


777=1 


M 

E ( M+1 ) 2 

777=1 

^ ✓ 

M(M+ 1) 2 


3 


Received Baseband Signal 


n = S bm + A 

with a 2 = S{\z\ 2 } = Nq. Again only the real part of the received 
signal is relevant for detection and we get 


tr — Sbm + 


with noise variance a 2 R = S{z^} = Nq/2. 

Decision Regions for ML Detection 


y/E~gd 

A 


m = 1 


mm 2 


- 4 -< 


m = M 


We observe that there are two different types of signal points: 

1. Outer Signal Points 

We refer to the signal points with rh = 1 and m = M as cmter 
signal points since they have only one neighboring signal point. 
In this case, we make on average 1/2 symbol errors if 

| Sbm\ A d\J Eg 

2. Inner Signal Points 

Signal points with 2 < m < M — 1 are referred to as inner 
signal points since they have two neighbors. Here, we make 
on average 1 symbol error if \tr — Sb m \ > dy/W g . 


Symbol Error Probability (SEP) 

The SEP can be calculated to 


Pm 


1 

M 
M - 1 


1 


(M - 2) + - • 2 


P (Vfl - Stm\ > dy/E^j 

M P (U Sbm ^ dy/Eg\ V [r R Sbm ^ d 


M- 1 
M 

M - 1 


( p ( 


r# - > d 


VX)+p{ 


T'R Sbm < d 


2P l r R - Sbm > dJ-p, 


= 2 


M 

M- 1 1 

M \J 71 Nq 


/ “ p (-5) * 




= 2 - 


.M - 1 1 
M y/2n 


I “ p ( _ b 


d2/ 


= 2^1 q 
M v 

Using the identity 


2rf2 l? 

iV 0 




y/X)) 


p m = 2 TzV I 6Es 


M 


(M 2 - l)JV 0 / ' 


we obtain 


We make the following observations 


1. For constant Es the error probability increases with increasing 
M. 

2. For a given SEP the required Es/Nq increases as 

101og 10 (M 2 — 1) « 201og 10 M. 

This means if we double the number of signal points, Le., M = 
2 k is increased to M = 2 fc+1 , the required Es/Nq increases 
(approximately) as 

20 log 10 (2 k+1 /2 k ) = 20 log 10 2 ss 6dB. 


Alternatively, we may express Pm as a function of the average 
energy per bit E\ )1 which is given by 

E = Es_ = E s 

h k log 2 M 

Therefore, the resulting expression for Pm is 


Pm = 2 - 


M-l 


M 


Q 


( / 61og 2 (M)E b \ 
[\l (M2 - 1 )Nq ) • 


An exact expression for the bit error probability (BEP) is more 
difficult to derive than the expression for the SEP. However, for 
high Es/Nq ratios most errors only involve neighboring signal 
points. Therefore, if we use Gray labeling we make approximately 
one bit error per symbol error. Since there are log 2 M bits per 
symbol, the PEP can be approximated by 


Pb 


1 

log 2 M 


Pm- 
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4.2.3 M - ary PSK 


■ For 2PSK the same SEP as for 2PAM results. 


P2 = Q 



■ For 4PSK the SEP is given by 


P 4 = 2Q 




■ For optimum detection of M - ary PSK the SEP can be tightly 
approximated as 


Pm 


2 Q 


( l 2\og 2 (M)E b . tt\ 

[\l No Sm M J • 


[ 
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■ The approximate SEP is illustrated below for several values of M. 

For M = 2 and M = 4 the exact SEP is shown. 



Eb/N 0 [dB] 


4.2.4 M- ary QAM 


■ For M = 4 the SEP of QAM is identical to that of PSK. 

■ In general, the SEP can be tightly upper bounded by 


Pm A 4Q 


( / 31og 2 (M)E h \ 
\\l (M - 1)N 0 ) ' 


■ The bound on SEP is shown below. For M = 4 the exact SEP is 
shown. 


188 | 



E h /N 0 [dB] 


4.2.5 Upper Bound for Arbitrary Linear Modulation Schemes 

Although exact (and complicated) expressions for the SEP and BEP of 
most regular linear modulation formats exist, it is sometimes more con- 
venient to employ simple bounds and approximation. In this sections, 
we derive the union upper bound valid for arbitrary signal constella- 
tions and a related approximation for the SEP. 

■ We consider an M - ary modulation scheme with M signal points 
s m: 1 < m < M, in the signal space. 

■ We denote the pairwise error probability of two signal points s ^ 
and s U: fi 7 ^ v by 

PEP(s Ai — » s u ) = P(s^ transmitted, s v , detected) 

■ The union bound for the SEP can be expressed as 


^ M M 

Pm ~ ~M 53 X/ PEP ( S M “ > 

/X=l v=l 


which is an upper bound since some regions of the signal space 
may be included in multiple PEPs. 

The advantage of the union bound is that the PEP can be usu- 
ally easily obtained. Assuming Gaussian noise and an Euclidean 
distance of d^ v = ||s ;/ — s u \\ between signal points s ^ and s v , we 
obtain 

mm / 

v 

Assuming that each signal point has on average Cm nearest neigh- 
bor signal points with minimum distance d m - m = min^^d^} and 
exploiting the fact that for high SNR the minimum distance terms 
will dominate the union bound, we obtain the approximation 



Note: The SEP approximation given above for MPSK can be 
obtained with this approximation. 

For Gray labeling, approximations for the BEP are obtained from 
the above equations with P& « Pm/ log 2 (M) 



[ 
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4.2.6 Comparison of Different Linear Modulations 

We compare PAM, PSK, and QAM for different M. 


m M = 4 
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uM = 16 



E h /N 0 [dB] 


■ M = 64 



E b /N 0 [dB] 


Schober: Signal Detection and Estimation 
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■ Obviously, as M increases PAM and PSK become less favorable 
and the gap to QAM increases. The reason for this behavior is the 
smaller minimum Euclidean distance d m - m of PAM and PSK. For 
a given transmit energy d m - m of PAM and PSK is smaller since the 
signal points are confined to a line and a circle, respectively. For 
QAM on the other hand, the signal points are on a rectangular 
grid, which guarantees a comparatively large d m - m . 

4.3 Receivers for Signals with Random Phase in AWGN 
4.3.1 Channel Model 

■ Passband Signal 

We assume that the received passband signal can be modeled as 


where z(t) is complex AWGN with power spectral density <& zz (f) = 
No, and </> is an unknown, random but constant phase. 

— (j) may originate from the local oscillator or the transmission 
channel. 

— <j) is often modeled as uniformly distributed, he., 


r(t) = V2Re { (e 3(t> s hm (t) + z(t)) e j2nfct } , 



[ 
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1/(2tt) 







— 7T 

7T (f 


■ Baseband Signal 

The received baseband signal is given by 

n(t) = e j4> s hm (t) + z(t). 

■ Optimal Demodulation 

Since the unknown phase results just in the multiplication of the 
transmitted baseband waveform Sb m {t) by a constant factor e J< ^, 
demodulators that are optimum for 0 = 0 are also optimum for 
<f> ^ 0. Therefore, both correlation demodulation and matched- 
filter demodulation are also optimum if the channel phase is un- 
known. The demodulated signal in interval kT < t < (k + 1)T 
can be written as 

r b [k\ = e ] ' t, s bm [k} + z[k\, 

where the components of the AWGN vector z[k\ are mutually 
independent, zero-mean complex Gaussian processes with variance 
of = N 0 . 


z(t) 


Sbm(t) 


n{t ) . 


-\ r b[k] r— 
ail ►Det 


m[k] 


[ 
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4.3.2 Noncoherent Detectors 

In general, we distinguish between coherent and noncoherent detec- 
tion. 

1. Coherent Detection 

■ Coherent detectors first estimate the unknown phase <f using 
e.g. known pilot symbols introduced into the transmitted signal 
stream. 

■ For detection it is assumed that the phase estimate <f is perfect , 
i.e., </> = </>, and the same detectors as for the pure AWGN 
channel can be used. 



z[k] 



■ The performance of ideal coherent detection with f = <f con- 
stitutes an upper bound for any realizable non-ideal coherent 
or noncoherent detection scheme. 

■ Disadvantages 

— In practice, ideal coherent detection is not possible and the 
ad hoc separation of phase estimation and detection is sub- 
optimum. 
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— Phase estimation is often complicated and may require pilot 
symbols. 

2. Noncoherent Detection 

■ In noncoherent detectors no attempt is made to explicitly es- 
timate the phase (f>. 
u Advantages 

— For many modulation schemes simple noncoherent receiver 
structures exist. 

— More complex optimal noncoherent receivers can be de- 
rived. 

4.3.2. 1 A Simple Noncoherent Detector for PSK with Dif- 
ferential Encoding (DPSK) 

■ As an example for a simple, suboptimum noncoherent detector we 
derive the so-called differential detector for DPSK. 

■ Transmit Signal 

The transmitted complex baseband waveform in the interval kT < 
t < (k + 1)T is given by 

Sbm(t) = b[k\g(t - kT), 

where g(t) denotes the transmit pulse of length T and b [k] is the 
transmitted PSK signal which is given by 

b[k] = e jQ M. 

The PSK symbols are generated from the differential PSK (DPSK) 
symbols a[k] as 


b[k] = a[k]b[k — 1], 


[ 
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where a[k\ is given by 

a[k] = e? Ae M, A©[fc] = 2 tt ( m - 1 )/M, m e {1, 2, . . . , M}. 

For simplicity, we have dropped the symbol index m in A0[fc] 
and a[k\. Note that the absolute phase ©[&] is related to the 
differential phase A0[/c] by 

Q[k\ = Q[k — 1] + A0[fc]. 


■ For demodulation, we use the basis function 
f(t) = —7¥r9{t) 

v-N 


■ The demodulated signal in the Arth interval can be represented as 


r' b [k] = e?* s/Egb[k] + z'[k\, 


where z'[k\ is an AWGN process with variance cr 2 ,, = No. It is 
convenient to define the new signal 


n[k] 



e ^ b[k] + z[k\ : 


where z[k\ has variance a 2 = Nq/E q . 



Differential Detection (DD) 

— The received signal in the kth and the (k— l)st symbol intervals 
are given by 

r b [k\ = e^a[k\b[k — 1] + z[k\ 

and 

r b [k — 1] = e^b[k — 1] + z[k — 1], 

respectively. 

— If we assume z[k\ « 0 and z[k — 1] ~ 0, the variable 

d[k] = r b [k]r* b [k - 1] 

can be simplified to 

d[k] = r b [k]r* b [k- 1] 

« e^a[k]b[k — 1 \{e^a[k\b[k — 1])* 

= a[k]\b[k]\ 2 
= a[k\. 

This means d[k\ is independent of the phase and is suitable 
for detecting a[k\. 

— A detector based on d[k\ is referred to as (conventional) dif- 
ferential detector. The resulting decision rule is 

m[k\ = argmin {| d[k\ — a^[A:]| 2 } , 

fh 

where a^[/c] = e- ?27r ( m_1 )/ M , 1 < fh < M. Alternatively, we 
may base our decision on the location of d[k\ in the signal space, 
as usual. 


d[k] 


m[k ] 


r b [k\ 



f 


Comparison with Coherent Detection 

- Coherent Detection (CD) 

An upper bound on the performance of DD can be obtained 
by CD with ideal knowledge of (j>. In that case, we can use 
the decision variable r&[&] and directly make a decision on 
the absolute phase symbols b[k\ = e J ’ e ^ £ {e = 
1, 2, . . . , M}. 

b[k] = argmin j \b[k] — e _J< ^rft[A:]| 2 | , 

m 1 J 

and obtain an estimate for the differential symbol from 

arh[k\ = b[k\ • b*[k - 1 ]* 

b[k] and b[k\ = £ | e j27r(m-i)/M| m = 1, 2, . . . , M} de- 

note the estimated transmitted symbol and a trial symbol, re- 
spectively. The above decision rule is identical to that for PSK 
except for the inversion of the differential encoding operation. 


e~M 



Since two successive absolute symbols b[k\ are necessary for 
estimation of the differential symbol a^k] (or equivalently ra), 
isolated single errors in the absolute phase symbols will lead to 
two symbol errors in the differential phase symbols, he., if all 
absolute phase symbols but that at time ko are correct, then all 
differential phase symbols but the ones at times ko and ko + l 
will be correct. Since at high SNRs single errors in the absolute 
phase symbols dominate, the SEP SEPj=jp SK of DPSK with CD 
is approximately by a factor of two higher than that of PSK 
with CD, which we refer to as SEPpg K . 

oppCD ^ ocupCD 
^ r DPSK ~ z ^ r PSK 

Note that at high SNRs this factor of two difference in SEP 
corresponds to a negligible difference in required E^/Nq ratio 
to achieve a certain BER, since the SEP decays approximately 
exponentially. 

- Differential Detection (DD) 

The decision variable d[k\ can be rewritten as 

d[k\ = r h [k]rl[k - 1] 

= (e^a[k]b[k — 1] + z[k])(e^b[k — 1] + z[k — 1])* 

= a[k\ + e~^b*[k — 1 \z[k\ + e^b[k\z*[k — 1] + z[k]z*[k — 1| , 

z eS ] 

where z e s[k\ denotes the effective noise in the decision vari- 
able d[k\. It can be shown that z e s[k] is a white process with 


variance 


2 _ 0 2 , 4 

'Vff - 2 < 7 z+^ 


= 2— + f — 






For high SNRs we can approximate <j\ as 


« 2az 
9 No 
e: 


- Comparison 

We observe that the variance cr? of the effective noise in the 
decision variable for DD is twice as high as that for CD. How- 
ever, for small M the distribution of z e s[k\ is significantly dif- 
ferent from a Gaussian distribution. Therefore, for small M 
a direct comparison of DD and CD is difficult and requires a 
detailed BEP or SEP analysis. On the other hand, for large 
M > 8 the distribution of z e s[k] can be well approximated 
as Gaussian. Therefore, we expect that at high SNRs DPSK 
with DD requires approximately a 3 dB higher E^/Nq ratio to 
achieve the same BEP as CD. For M = 2 and M = 4 this dif- 
ference is smaller. At a BEP of 10 -5 the loss of DD compared 
to CD is only about 0.8 dB and 2.3 dB for M = 2 and M = 4, 
respectively. 



.3.2.2 Optimum Noncoherent Detection 

■ The demodulated received baseband signal is given by 

r h = e J(t) s hm + 2 . 

■ We assume that all possible symbols are transmitted with 
equal probability, he., ML detection is optimum. 

■ We already know that the ML decision rule is given by 

rri = argmax {p(r \ Sbjh) } . 

m 

Thus, we only have to find an analytic expression for p(r |s^). The 
problem we encounter here is that, in contrast to coherent detec- 
tion, there is an additional random variable, namely the unknown 
phase </), involved. 


[ 
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i We know the pdf of r b conditioned on both 0 and s b m- It is given 

by 

I In - Sfm|| 2N 


p(n\sbm, <t>) = 


1 

W) 


f exp - 


N n 


I We obtain p(r h \si„j,) from p{ri\s hni , <j>) as 

OO 

p(r b \s b fh) = J p(r b \s bfn , 0)p$(>)d0. 

— OO 

Since we assume for the distribution of 4> 


P*(0) = { ft’ 


—7 r < (f) < TV 

otherwise 


we get 


7T 

P(nl^6m) — J P(Tb\strmi 4») 


exp - 


I In - e J b gfcmll 2 

N„ 


VoU 2 tt J 

—7 r 

1 1 f ( |n || 2 + llstmll 2 - 2Re{n«e^S6A}\ 

A4 J" 


Now, we use the relations 

Re{n • e^syf,} 


Re{|n»S6m|e J V^} 

In • s bm\ cos {<t> + a), 


where a is the argument of r& • s With this simplihcation 
can be rewritten as 

iwi 2 +n^ir 


p(n\Sbrh) = 


(7rw 0 ; 


F exp 


Nn 


— J exp (J^\ r b • s&m|cos(</> + d</> 


Note that the above integral is independent of a since a is inde- 
pendent of <j) and we integrate over an entire period of the cosine 
function. Therefore, using the definition 


h(x) = — 


i / exp(2 


we finally obtain 


IMMkhnll 2 

N 0 


p{n\ s m) = (^pexp (- 


I() y jy I'f'b • Sbfh\ 


Note that Io(^) is the zeroth order modified Bessel function of 
the first kind. 



x ► 


■ ML Detection 

With the above expression for the noncoherent ML de- 

cision rule becomes 


m = argmax {p(r \ Sbm ) } . 

rh 

= argmax {ln[p(r|s;,*)]}. 

fh 

lkw.ll 2 


= argmax 


N 0 


+ In 


Io ' 1% \ r b 9S bm\ 


205 | 

■ Simplification 

The above ML decision rule depends on No, which may not be 
desirable in practice, since No (or equivalently the SNR) has to be 
estimated. However, for x 1 the approximation 

ln[I 0 (x)] « x 

holds. Therefore, at high SNR (or small iVo) the above ML metric 
can be simplihed to 

m = argmax { — ||s^|| 2 + 2|r& • s^|} , 

rh 

which is independent of Nq. In practice, the above simplihcation 
has a negligible impact on performance even for small arguments 
(corresponding to low SNRs) of the Bessel function. 

4. 3. 2. 3 Optimum Noncoherent Detection of Binary Orthog- 

onal Modulation 

■ Transmitted Waveform 

We assume that the signal space representation of the complex 
baseband transmit waveforms is 

Sbi = [y/E, 0] T 

Sb2 = [0 VEbf 

■ Received Signal 

The corresponding demodulated received signal is 
r b = [<Py/E b + Zl z 2 ] T 


and 


n = [zi 


e^^/F b + z 2 f 
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if s b i and s b2 have been transmitted, respectively. z\ and z 2 are 
mutually independent complex Gaussian noise processes with iden- 
tical variances cr? = Nq. 

■ ML Detection 

For the simplification of the general ML decision rule, we exploit 
the fact that for binary orthogonal modulation the relation 

IMP = IMP- n 


holds. The ML decision rule can be simplified as 


m = argmax 


ikw.11 2 

N 0 


> Eb 1 

argmax < — — — + In 

Nq 


+ In 


Io ( N^ rb * Sbffl \ 


1() ( jy I'f'b • Sbm 


argmax < In 

m f 

argmax {\r b • s bfh \} 

m 

argmax {|n • M 2 } 


l0 ( \ r b • Sfoiil 


We decide in favor of that signal point which has a larger correla- 
tion with the received signal. 

■ We decide for m = 1 if 

\r b 9 s b i\ > \r b 9S b2 \. 


Using the definition 


r b = [r b i r b2 f, 


[ 


we obtain 


\r b { 


~y/W 

2 

~ni 


0 1 

\j h 2_ 

• 

0 

> 


• 

.%/aJ 



In 

2 > 

n| 2 . 




In other words, the ML decision rule can be simplihed to 

rh = 1 if |ri| 2 > |r2 1 2 

m = 2 if |ri| 2 < |r 2 | 2 
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Example: 

Binary FSK 

— In this case, in the interval 0 < t < T we have 


- If s bl (t) and s h2 (t) 
/iW = 

h{t) = 


Sbl{t) = 

s b2 (t) = 

are orthogonal, the basis functions are 
7r« 

0, otherwise 

f _J_pj27rA/t 

J VT e 

l 0 - 



0 <t<T 
otherwise 


[ 


208 | 


Receiver Structure 



t = T 


If sia(t) was transmitted, we get for r&i 

T 

m = / n(t)ft(t)dt 


0 

T 

J e?*sia(t)fi(t) 

0 

T 

^ 7 fi Sb2{t) 

o 

T 

J e^^dt + z, 


dt + z\ 


dt + z\ 


rh = 1 
rh = 2 


= \fhd.( J0 1 gj2xA/f , 

V 6 j2ttA/T + 1 

= ^ eJ0 j27rA/r ® J7rA/T ( eJ,rA/T " e_J7rA/T ) + 
= sinc(7rA/T) + ^ 




For the proposed detector to be optimum we require orthogo- 
nality, La, r&i = z\ has to be valid. Therefore, 

sinc(7rA fT) = 0 


is necessary, which implies 

A/r = i + ^, k €{..., -1, 0, 1, . . .}. 


This means that for noncoherent detection of FSK signals the 
minimum required frequency separation for orthogonality is 


A/ = 


1 

T' 


Recall that the transmitted passband signals are orthogonal 
for 


A/ = 


1 

2T’ 


which is also the minimum required separation for coherent 
detection of FSK. 


— The BEP (or SEP) for binary FSK with noncoherent detection 
can be calculated to 




[ 


210 | 



4. 3. 2. 4 Optimum Noncoherent Detection of On-Off Keying 

■ On-off keying (OOK) is a binary modulation format, and the 
transmitted signal points in complex baseband representation are 
given by 

Sbi = V^Eb 
Sb2 = 0 

■ The demodulated baseband signal is given by 

_ f e^y/2 Eb + z if m = 1 

rb \ z if m = 2 

where z is complex Gaussian noise with variance = Nq. 


[ 


■ ML Detection 

The ML decision rule is 

m = argmax < 

m l 


kU 2 

N 0 


+ In 


Io 
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■ We decide for m = 1 if 


2£* 

'No 


+ ln 


In 


f 0 

f 0 (-y/mn l) 


\/“L ( O, 


> ln(lo(0)) 

2^6 

N 0 


\n\ > 



Therefore, the (approximate) ML decision rule can be simplified 
to 


m = 1 if |r;,| > yT 

m = 2 if |n| < y^ 


■ Signal Space Decision Regions 


\: 171=1 

At 4 • 

TO = 2 / ■ 


.3.2.5 Multiple-Symbol Differential Detection (MSDD) of 
DPSK 

■ Noncoherent detection of PSK is not possible, since for PSK the 
information is represented in the absolute phase of the transmitted 
signal. 

■ In order to enable noncoherent detection differential encoding is 
necessary. The resulting scheme is referred to as differential PSK 
(DPSK). 


The (normalized) demodulated DPSK signal in symbol interval k 
is given by 

r b [k\ = e^b[k] + z[k] 

with b[k\ = a[k]b[k — 1] and noise variance a 2 = Nq/Es 7 where 
Es denotes the received energy per DPSK symbol. 

Since the differential encoding operation introduces memory, it is 
advantageous to detect the transmitted symbols in a block-wise 
manner. 

In the following, we consider vectors (blocks) of length TV 
r b = e^b + z, 

where 


n = [n[k\ r b [k - 1] . 

h 

1 

1 

,»sa. 

if 

b = [b[k\ b[k — 1] . . 

1 

1 

z = [z[k\ z[k — 1] . . 

1 

1 

O 


Since z is a Gaussian noise vector, whose components are mutu- 
ally independent and have variances of a 2 = TV 0 , respectively, the 
general optimum noncoherent ML decision rule can also be applied 
in this case. 

In the following, we interpret b as the transmitted baseband signal 
points, i.e., s b = 6 , and introduce the vector of TV — 1 differential 
information bearing symbols 

a = [a[k\ a[k — 1] ... a[k — (TV — 2)]] T . 

Note that the vector b of absolute phase symbols corresponds to 
just the TV — 1 differential symbol contained in a. 


The estimate for a is denoted by a, and we also introduce the 
trial vectors a and b. With these definitions the ML decision rule 
becomes 


a = argmax 


= argmax < — — — + In 
Nq 


mi 

N 0 

N 


+ In 




N 0 


|n»&| 


argmax {ln«b|} 

N-l 

rb[k — v]b*[k — v\ 


= argmax 


In order to further simplify this result we make use of the relation 

N-2 

b[k — v\ = Y[ — l$\k — (N — 1)] 

(JL— V 

and obtain the final ML decision rule 


a = argmax 


= argmax 


N - 1 


N-2 


yt n[k - H LI “ m] N[k - (N - 1)] 

w—0 n=v 

N-l N-2 

yi r ^ k - h n ~ 


Note that this hnal result is independent of b[k — (N — 1)]. 

We make a joint decision on N — 1 differential symbols a[k — v\ 
0 < v < N — 2, based on the observation of N received signal 
points. N is also referred to as the observation window size. 


The above ML decision rule is known as multiple symbol differ- 
ential detection (MSDD) and was reported first by Divsalar & 
Simon in 1990. 


In general, the performance of MSDD increases with increasing 
N. For N — ■» oo the performance of ideal coherent detection is 
approached. 

N = 2 

For the special case N = 2, we obtain 


= argmax {|r&[fc]| 2 + | r b [k - 1]| 2 + 2Re{r b [k]a*[k]rl[k - 1]}} 


= argmin {\n[k]r b [k - 1]| 2 + |a[£;]| 2 - 2Re{r b [k]a*[k]rl[k - 1]}} 


Obviously, for N = 2 the ML (MSDD) decision rule is identical to 
the heuristically derived conventional differential detection decision 
rule. For N > 2 the gap to coherent detection becomes smaller. 


a 


\k\ 


argmax 

a[k\ 



argmax {| r b [k]a*[k] + r b [k - 1] | } 

a[k] 



a[k] 

argmin {— 2Re{r b [k]a*[k\rf[k — 1]}} 

a[k] 


a[k ] 

argmin {| d[k\ — a[A:]| 2 } 

a[k] 


with 


d[k\ = r b [k]rl[k - 1]. 


[ 


■ Disadvantage of MSDD 

The trial vector a has N — 1 elements and each element has M 
possible values. Therefore, there are M N ~ X different trial vectors 
a. This means for MSDD we have to make M n ~ 1 /(N — 1) tests 
per (scalar) symbol decision. This means the complexity of MSDD 
grows exponentially with N. For example, for M = 4 we have to 
make 4 and 8192 tests for N = 2 and N = 9, respectively. 

■ Alternatives 

In order to avoid the exponential complexity of MSDD there are 
two different approaches: 

1. The first approach is to replace the above brute-force search 
with a smarter approach. In this smarter approach the trial 
vectors a are sorted first. The resulting fast decoding algorithm 
still performs optimum MSDD but has only a complexity of 
IVlog(fV) (Mackenthun, 1994). For fading channels a similar 
approach based on sphere decoding exists (Pauli et ah). 

2. The second approach is suboptimum but achieves a similar 
performance as MSDD. Complexity is reduced by using de- 
cision feedback of previously decided symbols. The resulting 
scheme is referred to as decision-feedback differential detection 
(DF-DD). The complexity of this approach is only linear in N 
(Edbauer, 1992). 


BER 


■ Comparison of BEP (BER) of MSDD (MSD) and DF-DD for 
4DPSK 



10log 10 (E b /N 0 ) [dB] 


Schober: Signal Detection and Estimation 
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4.4 Optimum Coherent Detection of Continuous Phase Mod- 
ulation (CPM) 

■ CPM Modulation 

— CPM Transmit Signal 

Recall that the CPM transmit signal in complex baseband rep- 
resentation is given by 


where I is the sequence {/[&]} of information bearing symbols 
I[k\ G {±1, ±3, . . . , ±(M — 1)}, (j)(t, I) is the information 
carrying phase, and <f o is the initial carrier phase. Without 
loss of generality, we assume </>o = 0 in the following. 

- Information Carrying Phase 

In the interval kT < t < (k + 1 )T the phase <f>(t, I) can be 
written as 


where Q[k\ represents the accumulated phase up to time kT, h 
is the so-called modulation index, and q(t) is the phase shaping 
pulse with 



L—l 


<j)(t, I) = Q[k] + 27 rh I[k — u\q(t — [k — v]T) 


+27 ihl[k\q(t — kT), 



0 , 

monotonic, 

1 / 2 , 


t < 0 

0 < t < LT 
t > LT 


— Trellis Diagram 

If h = q/p is a rational number with relative prime integers 
q and p, CPM can be described by a trellis diagram whose 
number of states S is given by 


with complex AWGN z(t), which has a power spectral density of 
*zz(f) = N 0 . 

ML Detection 

— Since the CPM signal has memory, ideally we have to observe 
the entire received signal r&(£), — oo < t < oo, in order to 
make a decision on any I[k] in the sequence of transmitted 
signals. 

— The conditional pdf p(r&(t)|s&(£, I)) is given by 


where / E {±1, ±3, . . . , ±(M — 1)} is a trial sequence. For 



Received Signal 

The received complex baseband signal is given by 


n(t) = s b (t, I) + z(t ) 


p(n(t)\sb(t, i)) oc exp 



ML detection we have the decision rule 


I = argmax {p(r b (t)\s b (t, /))} 

= argmax {ln\p(r b (t)\s b (t, I))]} 

i 

= argmax | - J \ r b (t) - s b (t, J)| 2 d< j 

! oo oo 

- J \n(t)\ 2 dt - J | s b (t, 

— oo — oo 

oo ^ 

+ J 2Re{r b (t)s b (t, /)} dt 

— OO J 

= argmax II Re{r b (t)s b (t, J)} dt, 1 


ni 2 df 


where 7 refers to the ML decision. If 7 is a sequence of length 
K , there are M K different sequences 7. Since we have to cal- 
culate the function J ^ R e{r&(£)sj(£, 7)}d£ for each of these 
sequences, the complexity of ML detection with brute-force 
search grows exponentially with the sequence length iL, which 
is prohibitive for a practical implementation. 


Viterbi Algorithm 

— The exponential complexity of brute-force search can be avoided 
using the so-called Viterbi algorithm (VA). 

— Introducing the definition 

(k+l)T 

A[fc] = J Re{r b (t)sl(t, I)}dt, 

— OO 

we observe that the function to be maximized for ML detection 
is A[oo]. On the other hand, A [k] may be calculated recursively 
as 

kT (k+l)T 

A[fc] = J Re{r b (t)s* b (t, I)}dt + J Re{r b (t)s* b (t, I)} dt 


(k+l)T 


A[fe] = A[fe — 1] + J R e{r b (t)s* b (t, I)} 

kT 

A [k\ = A [k — 1] + A [k\, 
where we use the definition 

(k+l)T 

A[fc] = J Re{r b (t)s* b (t, I)}dt. 


d t 


(k+l)T 


J 


Re<! rb(t) exp -j 


0[fc] + 27 ih I[k — v\ 


• q(t — [k — v\T) + 27 rhl[k\q(t — kT) 


d t. 


A [k\ and A [A;] are referred to as the accumulated metric and 
the branch metric of the VA, respectively. 

— Since CPM can be described by a trellis with a finite number 
of states S = pM L ~ x ( S = 2 at time kT we have to 
consider only S different A[S[fc — 1]. k — 1 ], Each A[£[£; — 
1], k — 1] corresponds to exactly one state S[k — 1] which is 
defined by 

S[k - 1] = [Q[k - 1], I[k — (L — 1)], . . . , I[k - 1]]. 


S[k - 2] S[k - 1] S[k\ 

1 • • • 









For the interval kT < t < (k + 1)T we have to calculate 
M branch metrics X[S[k — 1], I[k], k] for each state S[k — 1] 
corresponding to the M different I[k\. Then at time ( k + 1)T 
the new states are defined by 

S[k] = [6[fc], I[k - (L - 2)], ... , I[k}\, 

with Q[k] = Q[k — 1] +7vhl[k — (L — 1)]. M branches defined 
by S[k — 1] and I[k] emanate in each state S[k]. We calculate 


the M accumulated metrics 

A [S[k - 1], I[k], k] = A[S[/c - 1], k - 1] + A [S[k - 1], I[k], k] 
for each state S[k\. 



From all the paths (partial sequences) that emanate in a state, 
we have to retain only that one with the largest accumulated 
metric denoted by 

A[S[fc], k] = max /[|fc| {A [S[k - 1], I[k], A:]} , 

since any other path with a smaller accumulated metric at time 
(k + l)T cannot have a larger metric A[oo]. Therefore, at time 
(k + 1)T there will be again only S so-called surviving paths 
with corresponding accumulated metrics. 



— The above steps are carried out for all symbol intervals and at 
the end of the transmission, a decision is made on the trans- 
mitted sequence. Alternatively, at time kT we may use the 
symbol I[k — ko] corresponding to the surviving path with the 
largest accumulated metric as estimate for I[k — ko]. It has 
been shown that as long as the decision delay ko is large enough 
(e.g. ko > 5 log 2 (A)), this method yields the same performance 
as true ML detection. 

- The complexity of the VA is exponential in the number of 
states, but only linear in the length of the transmitted se- 
quence. 

■ Remarks: 

— At the expense of a certain loss in performance the complexity 
of the VA can be further reduced by reducing the number of 
states. 

— Alternative implementations receivers for CPM are based on 
Laurent’s decomposition of the CPM signal into a sum of PAM 
signals 


